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We consider a driven damped anharmonic oscillator which 
classically leads to a bistable steady state and to hysteresis. 
The quantum counterpart for this system has an exact an- 
alytical solution in the steady state which does not display 
any bistability or hysteresis. We use quantum state diffusion 
theory to describe this system and to provide a new perspec- 
tive on the lack of hysteresis in the quantum regime so as to 
study in detail the quantum to classical transition. The anal- 
ysis is also relevant to measurements of a single periodically 
driven electron in a Penning trap where hysteresis has been 
observed. 

PACS 03.65.Bz, 42.50.Lc, 05.30Ch 



I. INTRODUCTION 

Emergence of classical properties through interaction 
with the environment has been the subject of extensive 
studies. In this context, a study of simple open quantum 
systems can provide a clue to understanding the mech- 
anism of the quantum-classical transition. In this paper 
we wish to study one of the simplest nonlinear quantum 
systems - the anharmonic oscillator. This system, when 
it is damped and driven, exhibits different behaviour in 
the steady state when described either quantum mechani- 
cally or classically with the later showing bistability E10 . 
In order to describe the emergence of classicality, such 
differences have to be clarified and the two results recon- 
ciled. 

The anharmonic oscillator is of particular interest for 
several reasons. First, its simplicity allows any complex- 
ity related to the model to be avoided. Second, it is 
the archetypical model for dealing with nonlinearities in 
quantum mechanics and has been widely used to describe 
a great variety of systems. In particular, it was intro- 
duced by Drummond and Walls |3|] to describe dispersive 
optical bistability and more recently has been used by 
Bortman and Ron to study the relativistic motion 
of a resonantly driven electron in a Penning trap. 
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Drummond and Walls || derived the exact steady 
state expectations values of the photon distribution func- 
tion hence showing that in the quantum regime this 
model does not exhibit bistability or hysteresis. One 
would expect the classical or quantum results to be recov- 
ered in the appropriate limit. For instance, the level of 
excitation, (a^a), has been used to define such a limit M. 
For low excitation number the quantum result should ap- 
ply while for high excitation number the classical result 
is expected to give the correct behavior. This description 
is reasonable but it does not demonstrate how the transi- 
tion from classical to quantum occurs when the excitation 
number takes an intermediate value and whether the sys- 
tem suddenly becomes bistable or whether the bistability 
appears smoothly. 

Drummond and Walls ||] state that the extent to which 
bistability is observed will depend on the quantum fluc- 
tuations, which in turn determine the time for random 
switching between the two stable states. In this case the 
classical-quantum transition is of a statistical nature and 
the classical result appears in the limit where the switch- 
ing time can be considered very large compared to an 
observational time. 

Bortman and Ron presented a quantum-mecha- 
nical description of a single atomic system which is ac- 
cessible to experiment ||. They also state that if one 
does not deal with the effect of fluctuations upon the 
two stable states, that is, its influence upon the time scale 
of the stability, bistability is not destroyed by quantum 
fluctuations and should be observed even in the case of a 
low level of excitation. According to this description the 
system is bistable in both regimes and the observation of 
bistability thus depends on the experimental set up, with 
no fundamental restriction on its observation. 

Motivated by these controversial points of view of 
Drummond and Walls || and Bortman and Ron |l],^| in 
this paper the quantum-classical transition of the driven 
damped anharmonic oscillator is investigated with the 
help of the quantum state diffusion (QSD) method. In 
the context of typical quantum optical problems, for ex- 
ample, the QSD method describes the continuous moni- 
toring of the state of a photon source by individual pho- 
toelectric detection processes which involve heterodyning 
with a classical intense photon source However, the 
QSD method has also been proposed as a phenomenolog- 
ical theoretical description of arbitrary individual quan- 
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turn measurement processes [PHlTI- It has already been 
demonstrated that the QSD method, considered as a dy- 
namical theory for single quantum systems, is a valuable 
tool in the understanding of the emergence of classical 
chaos in open quantum systems | Q | . Thus it is expected 
that it will also be useful in obtaining new insights into 
the connection between classical and quantal behaviour 
of the driven damped anharmonic oscillator coupled to 
a reservoir, a physical system whose classical dynamics 
does not exhibit chaos. This is of particular interest 
here as this quantum system does not exhibit hystere- 
sis whereas the corresponding classical system does. It 
will be shown that QSD provides a mechanics for bistable 
motion in phase space which is consistent with the quan- 
tum mechanical steady state result that bistability does 
not appear for mean values over an ensemble. Further- 
more, QSD helps to understand the crucial role played by 
the physical time scales which characterize the approach 
to the classical equilibrium positions in phase space and 
the jumps between these two classically stable equilib- 
rium positions due to quantum fluctuations. In addition, 
the question is addressed as to what can be measured, 
in principle, in an experiment Q on such a system and 
under what conditions such an experimental observation 
will yield bistability and hysteresis. Finally the assertions 
of Drummond- Walls and Bortman-Ron are questioned in 
the light of QSD. 

The paper is organized as follows: In section [n|, the 
QSD model is briefly depicted. In section III the anhar- 
monic oscillator model is presented, and a brief review 
of the classical and quantum results related to bistability 
in the steady state is provided. Section IV describes the 
anharmonic oscillator from the point of view of QSD. Fi- 
nally, section ^ presents a discussion of the results and 
our conclusions. 



II. QUANTUM STATE DIFFUSION 

Open quantum systems are represented by the density 
operator p which evolves in time according to a mas- 
ter equation. The most general master equation in the 
Markov approximation, which preserves trace and pos- 
itivity of the density operator p can be written in the 
Lindblad form [fL3| 



\dip) = --H\ip)dt 



-[H,p]+J2(L j pL 



1 



1 



p l]lA (i) 



where H is the Hamiltonian and Lj are Lindblad oper- 
ators which represent the effects of the environment on 
the system. 

The quantum state diffusion model (QSD) represents 
one of the several possible unravellings of the master 
equation. According to the QSD model ]^- pr| , ^5| , ^6[ , 
open quantum systems are represented by pure states 
\ip), which describe individual systems. Evolution of the 
state \ip) is given by a Langevin-Ito differential equation 



j 

+ J2(L j -(L j )Mdt j . (2) 



The d£j are random differential variables representing in- 
dependent complex Wiener processes. They satisfy the 
following mean relationships 

AT(d&) = M(dfcd£fc) = 0, M(d&dft) = 6 jk dt. (3) 

M represents a mean over an ensemble and (Lj) — 
(ip\Lj\ip) the quantum expectation of the operator Lj 
in the pure state \ip). 

The QSD trajectories are compatible with the master 
equation in the sense that the ensemble average of the 
projector \ip)(ip\ reproduces the density operator p: 



p = M(MM). 



(4) 



Thus, expectations values (A) p — Tr(pA) of an operator 
A can be computed as the ensemble mean of the quantum 
expectations values (A)^ of the pure state \tp). 



III. ANHARMONIC OSCILLATOR 

We consider a driven anharmonic oscillator coupled to 
a thermal bath |^,|l7|-|l9| , the temperature of which is set 
to zero (T = 0) as the prototype model showing bista- 
bility in the classical domain. The damping of this oscil- 
lator with rate K, is described by the Lindblad operator 
L = y/Jia. The Hamiltonian in a frame rotating with the 
frequency u> of the driving field reads 



H = fiAuja^a + hf3(a) + a) + hxi^a) 2 



(5) 



where a 



(^) 1/2 Q + i(- 



_U/2p an d a t are the an- 
nihilation and creation operators related to the position 
Q and momentum P of the oscillator (m is the mass of 
the particle) . Here the parameter Auj — ujq — u> measures 
the detuning between the eigenfrequency of the oscillator 
and the driving force. The parameters f3 and x charac- 
terize the amplitude of the driving force the strength of 
the anharmonicity. In the following only positive values 
of x w iU be considered. 

This Hamiltonian whose corresponding classical dy- 
namics is integrable is well known. It has been used to 
describe various physical phenomena such as dispersive 
optical bistability ^,0, driven tunneling (2TJ and hys- 
teresis in an atomic system ^JJj| within the framework 
of the rotating wave approximation. In the context of 
optical bistability, for example, the QSD equation with 
Hamiltonian (5) and L — ^JUa describes the continuous 
monitoring of the electromagnetic field of frequency ujq by 



2 



individual photon detection processes which involve het- 
erodyning with a classical intense photon source. With- 
out the adiabatic approximation (i.e. for a time depen- 
dent driving 0) this physical system has been discussed 
as a model which exhibits quantum chaos |l2| . It has 
also been used to describe more fundamental aspects like 
the effect of non linearities in master equations , 11011- 
Markovian approximations fis}| , and more recently in a 
study of localization processes An appealing fea- 

ture of this system is that exact quantum results exist 
for the correlations function jjj, the spectrum, and even 
the dynamics jl^,^3| in the absence of driving. 

The classical equivalent for this system exhibits hys- 
teresis in the steady state while the quantum system 
does not. In the next two sections well known results re- 
garding the classical and quantum systems are presented 
in order to make this presentation self consistent. 



The first condition expresses the fact that the detuning 
has to be oriented in the right direction in order to com- 
bine its effect with that of the anharmonicity. The second 
shows that detuning and anharmonicity must be large 
enough in order to compensate dissipation, and the third 
condition gives limiting values for the driving strength. 

Using the theory of linear stability j^,[l) , it can be ver- 
ified that when three solutions are present, one of them 
is always unstable. Thus the domain of parameters is 
divided into two regions, one showing bistability and the 
other purely mono-stable behaviour (see figure ^|). In 
short, the classical model exhibits a bistable steady state 
when the parameters satisfy the above conditions leading 
to hysteresis. 

B. The quantum limit 



A. The classical limit 



The classical equation of motion can be obtained |3J by 
factorizing the quantum correlation functions (afar) — * 
(at)(a) 2 

= -i{/3+ (ALU + x)a + 2 X a 2 a*} - i«a (6) 

where a is the mean field amplitude a = (a) in the (semi- 
) classical limit. 

In the steady state regime, the excitation number \a\ 2 
can be obtained by solving the following equation 



M 2 = 



B 2 



( K /2) 2 + (A W + x + 2 X M 2 r 



(7) 



Once the excitation number \a\ 2 is known the real and 
complex part of the mean amplitude a can be derived 
using 



Re(a) _ Acj + x + 2 X |a| 2 



Im{a) 



at 
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Expression (0) shows that the classical anharmonic oscil- 
lator can display one, two or three solutions depending 
on the choice of parameter values (see figure |l|) . 
Provided the following conditions are fulfilled 



x(Alu + x) < 



Au; + x 



k/2 



> 



Using the complex P representation, Drummond and 
Walls H have solved the master equation for the density 
operator. They obtained an analytical expression for the 
moments (a) n a m ) p in the steady state. Notice that in this 
section, the brackets {■■ -) p represent the ensemble mean 
of the quantum mechanical expectation values. Their 
result reads 



Y(c)Y(c* 



2 J Y(c + m)Y(c* +n) 
F(c + m, c* + n, z) 



F(c, c*, z) 



(10) 



where Y is the Gamma function and F = 0F2 the gener- 
alized Gauss hypergeometric series pJ : 



oF 2 (c, d,z) = ^2 



Y(c)Y(d) 



-** n\ Y(c + n)Y(d + n) 



(ii) 



The coefficients c and z depend on the physical parame- 
ters in the following way c = (Aw + x)/x ~ * K /(2x) and 
z = 2(P/ X ) 2 . 

The mean excitation number is of particular interest 
here, it is given by 



(a t o) p 



/? 2 



F{c+ l,c* + l,z) 



(Auj + x) 2 + (k/2) 2 F(c,c*,z) 



(12) 



27 X P 2 



(Alo + x) 3 



3k/2 
Alu + x 



-1 2 



< 



k/2 
Au + x 



(9) 

n 3 



Using the properties of the hypergeometric series 0F2 one 
can show that the quantum result is never bistable and 
thus does not show any hysteresis (figure 0) . 



then equation (0) has three solutions. Outside of this 
range, only one solution is expected. 



IV. QSD FOR THE ANHARMONIC OSCILLATOR 

In this section, the problem is tackled using the quan- 
tum state diffusion (QSD) model. According to QSD, 
the equation of motion for the mean field amplitude (a) 
is 
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d(a) = -i [{Alo + x )(a) + (3 + 2x(a t a 2 }] dt - ^(a)dt 



+V^((a 2 ) - (a) 2 )de 
+VK((a t a) - (a t )(a))d^*. 



(13) 



In this equation, the expectation values are taken in the 
pure state describing the evolution along a quantum 
trajectory. If one wants to describe the evolution of the 
mean value (over an ensemble), one can take the mean 
on both sides of the equation (|l|) to obtain 



d(a), 
dt 



-i [{Alo + x){a) P + 13 + 2 X (aV)J - ^(a) f 



(14) 



One can check easily that factorizing the quantum cor- 
relations in both of the two preceding equations leads to 
the classical equation of motion (^|). Thus, neglecting 
the quantum correlations corresponds to ignoring over- 
lapping effects of the wave packet and quantum fluctua- 
tions. 



A. Simulation of the dissipative dynamics 

The parameters are chosen in the bistable domain. The 
dynamical evolution of the QSD equation ( [l3"| ) is com- 
puted numerically using the moving basis or mixed rep- 
resentation simulation method (MQSD) The evo- 
lution is computed over a period of time of the order of 
1/k, the dissipative time. In this situation, the system 
evolves toward a different "stationary" state depending 
on the chosen initial state. The QSD evolution shows 
two different limit points where all the trajectories tend 
to go after some transient dissipative time. These two 
points are called equilibrium points. 

When the trajectory starts far away from the two equi- 
librium points, it approaches one of them, depending in 
which basin of attraction it starts in, rotating with a fre- 
quency given by the detuning Alu. If the wave packet is 
initially spread out, it tends to localize jll| to a coherent 
state during the dissipative transient. Here the state is 
said to be localized if its spread is much smaller than the 
distance between the two equilibrium points. Thus the 
quantum fluctuations have less and less effect. 

After the transient damped motion towards one of the 
equilibrium points, the wave packet remains localized. 
This strong localization property allows one to describe 
the quantum system in a quasi classical way. 

As a consequence of the anharmonicity, the system 
does not however preserve the coherent states to which it 
is driven by the dissipative terms. Hence, the quantum 
correlations never vanish, and the quantum fluctuations 
act on the wave packet whose center will fluctuate around 
it's equilibrium point. The dynamics are now dominated 
by quantum fluctuations. 

These fluctuations have a mean frequency and an am- 
plitude which depends on the position of the equilibrium 



point in phase space. The fluctuations are bigger for 
the equilibrium point situated further away from the ori- 
gin. This is an effect due to dissipation, represented by 
L = ti/kcl, which produces a dynamical behaviour like 
(a) ~ —n(a) clearly attracting the system towards the 
origin (a) — and not to some local minima. 

At this point the QSD trajectories can be roughly seen 
as classical trajectories subjected to noise. There are two 
distinct equilibrium points leading to bistable behaviour 
similar to the classical one. 



B. Recovery of the quantum result 

We know that the QSD model reproduces the quantum 
result in mean over an ensemble but if QSD shows bista- 
bility, how can the quantum result be recovered ? The 
answer is that the evolution described above is stable 
over a very long time compared to the dissipative time, 
but if one integrates over a longer time, one sees that any 
trajectory goes from the neighborhood of one equilibrium 
point to the other. This transition happens in mean after 
a time called the transition time or exit time, which can 
be much longer than the dissipative time. 

In order to observe the transition between the two 
equilibrium points, the parameters are chosen such that 
the maximal excitation number is set at an intermediate 
value between the classical and quantum limits. Also, the 
integration is now carried over a long time, typically 10 2 
to 10 times the dissipative time. The system is initially 
set in a coherent state centered far away from the two 
fixed points and evolved in time with QSD. Figure || rep- 
resents a trajectory in phase space which shows the dis- 
sipative part of the trajectory, followed by a long period 
of fluctuations around the attracting equilibrium point. 
This first part of the time evolution of the trajectory, i.e. 
its approach to equili brium, has been described in the 
previous section IV A. If the integration is carried on, 



the trajectory will suddenly jump to the neighbourhood 
of the other equilibrium point. (The "jump" described 
here is a diffusive process which allows the quantum tra- 
jectory to go from the basin of attraction of one equilib- 
rium point to the other.) The system will remain around 
the second equilibrium point during some time and then 
come back to the first point. The time spent around 
each equilibrium point is such that the quantum expec- 
tation value (a) p is recovered in mean. Due to quantum 
fluctuations induced by the coupling to the reservoir the 
equilibrium points become metastable. 

This transition, which occurs in a time shorter than the 
dissipative time, might be viewed as a tunneling process. 
Savage and Cheng [gl| have investigated whether bista- 
bility can be associated with quantum superpositions of 
states in either well. They have introduced a distinction 
between coherent and diffusive mechanisms for quantum 
tunneling. In our simulations, the wave packet initially 
localized in one well becomes delocalized when it crosses 
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the barrier, making the distinction between these two 
mechanisms of tunneling artificial (see figure Q). Once 
the barrier is crossed, the wave packet localizes again. 

To confirm the previous description, a mean over an 
ensemble is considered. Figure ^represents the time evo- 
lution of the mean position and mean momentum. The 
mean is computed over 100 trajectories. The system is 
initially placed in a coherent state centered at the classi- 
cal equilibrium point. This point is unstable with respect 
to the other equilibrium point. The mean position and 
momentum evolve, roughly as an exponential decay, to 
the quantum stationary values given by the exact quan- 
tum result. For this typical example, the quantum result 
is very close to one of the equilibrium points, because the 
transition time from the initial to the final point is much 
shorter than in the opposite direction. 

Hence, the mean result confirms that the initial equi- 
librium point is unstable compared to the other one as 
a consequence of the quantum fluctuations. This rel- 
ative instability explains why the quantum description 
does not show bistability. According to QSD even in the 
quantum regime, the system is bistable, but the bista- 
bility is hidden by the fluctuations which make the wave 
packet move from a fairly localized state in one well to a 
localized state in the other well. 



V. DISCUSSION 

We have used QSD to describe the driven damped an- 
harmonic oscillator in an intermediate regime between 
quantum and classical. It has been shown that states 
localize along a quantum trajectory and a transition be- 
tween the two equilibrium points of the system has been 
observed. The localization gives a quantitative justifica- 
tion for the classical analogy in which a localized particle 
moves in a double well potential. This analogy has often 
been used on a purely qualitative level without any fur- 
ther justification. The transition between the two equi- 
librium points allows one to recover the quantum result 
and reconciles quantum and classical descriptions. It is 
worth emphasizing the following aspects: 

A. Transition time 

The QSD model, by introducing quantum fluctuations 
in an explicit way, shows explicitly how classical bista- 
bility disappears. Furthermore it introduces a new time 
scale, the transition time characteristic of the transition 
between the two (classical) equilibrium points. More pre- 
cisely there are two transitions times, the transition time 
from one equilibrium point to the other and a different 
time associated with the reverse transition. Because one 
of these transition times is in general much smaller than 
the other one, the transition considered here starts from 



the less stable equilibrium point which is located further 
away from the origin in phase space. 

In a sequence of papers, Vogel and Riskcn j2^] have cal- 
culated the transition rates by solving the equations of 
motion for quasidistribution functions using the matrix 
continued-fraction method. They also obtained analyt- 
ical results for the transition rate in the limit of large 
excitation numbers and low damping. 

The rate of decay shown on the figure || is an approxi- 
mation to the mean transition rate between the two equi- 
librium points. The time needed for the decay is clearly 
much larger than the dissipative time. 

Thus QSD not only gives a qualitative description but 
can also be used easily to obtain numerical estimates of 
the relevant transition times. We will not address any 
further the question of the determination of the transition 
time in this paper as one can use the accurate results of 
Vogel and Risken which confirm the possibility of very 
large transition times compared to the dissipative time. 

B. Ideal experiment 

Let us consider the following ideal experiment : (see 
for a practical realization) In order to see hysteresis a 
single quantum mechanical anharmonic oscillator is mea- 
sured continuously under conditions in which its classical 
counterpart would be bistable. Let us assume that the 
driving frequency u) is varied step by step from low to 
high frequencies and reversed, spanning twice the classi- 
cally bistable domain. Once the frequency is modified, 
the experimenter waits a time, called the measurement 
delay t m , which is assumed to be much longer than the 
dissipative time, before doing any further measurement. 
The excitation number of the oscillator is measured be- 
fore changing the frequency of the driving force again. 
Thus this type of resonance experiment corresponds to 
an adiabatic sweeping of the frequency with respect to 
the "fast" dissipative dynamics. Furthermore, let us as- 
sume that this continuous measurement is ideal in the 
sense that the excitation number of the harmonic oscil- 
lator can be measured nondestructively and that it does 
not perturb the systems transition rates, i.e. it is a non- 
demolition measurement. Such a measurement can be 
performed, for example, in an optically bistable system 
by monitoring the state of the field mode by heterodyning 
with an intense classical photon source. Alternatively, 
this measurement can also be realized by observing the 
relativistic motion of a resonantly driven electron in a 
Penning trap, like in the recently performed experiment 
of Ref. @. 

According to QSD such a measurement should repro- 
duce the curve depicted in figure |^ showing hysteresis. 
The experimenter should obtain such a curve fluctuat- 
ing around one of the two classical steady state values 
for a while and then jumping to the other value. The 
combination of the two jumps occurring when the driv- 
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ing frequency is ramped from low to high frequencies and 
reversed allows one to define the detuning width Af2 as 
the size of the bistable region. Figure || represents such 
a result and shows the detuning width Af2 for this par- 
ticular realization. The detuning width is different for 
each realization of this experiment, the transition being 
a stochastic event. 

An experiment carried out in the classical limit does 
not show any fluctuations and the two jumps occur al- 
ways at the same detuning value. In this case the detun- 
ing width Afl corresponds to the full size of the bistable 
region. 

If one uses the density matrix to describe such an ideal 
experiment, the result will also show two distinct transi- 
tions. The mean detuning width depends not only on the 
characteristic physical parameters of the driven damped 
anharmonic oscillator but also on the measurement de- 
lay t m . If the mean transition time r is much larger than 
the dissipative time, i.e. t 3> one can distinguish 
between the two limiting cases: (1) If the measurement 
delay is small relative to the transition time, i.e. r ^> t rn , 
then the mean detuning width Af2 has a finite value, 
showing bistability. (2) At the opposite extreme, i.e. for 
tm 2> t, the detuning width is equal to zero, showing no 
hysteresis at all, in agreement with the quantum steady 
state result. 



C. Bistability and the classical limit 

The classical limit is valid for high excitation number. 
In this limit the mean transition time is so large and 
the transition between the two equilibrium points so in- 
frequent that they can be neglected. In this limit any 
finite observational time satisfies the conditions for the 
experimental observation of hysteresis. 

When the excitation number decreases sufficiently for 
the mean transition time to take a very large but accessi- 
ble value, one has to distinguish between the cases where 
the measurement delay is smaller or larger than the tran- 
sition time. The classical description is still valid in the 
former case but does not apply anymore in the latter. 
The fluctuations have to be taken into account for a cor- 
rect description of this situation. 

If we continue lowering the excitation number, still 
keeping the mean transition time large compared to the 
dissipative relaxation time, then the classical theory no 
longer gives a good description of the dynamics since even 
for a measurement delay much smaller than the transi- 
tion time it predicts a fixed detuning width. If one uses 
quantum theory, which includes the fluctuations, one will 
be able to obtain the correct behaviour. 

Finally when the excitation number is small, the clas- 
sical theory is no longer valid. One has to use the quan- 
tum theory and specify the measurement delay in order 
to describe correctly the result of an experiment. The 
bistability is not destroyed in any of these cases but it is 



simply hidden by the quantum fluctuations. 

This situation is very similar to that of a classical 
driven anharmonic oscillator coupled with a thermal bath 
with non-zero temperature. Introducing thermal fluctu- 
ations also hides the bistability of the steady state and 
introduces a classical transition time (see |25|1 ). In or- 
der to observe bistability, one has to introduce a mea- 
surement delay much larger than the relaxation time in 
the absence of thermal fluctuations, but shorter than the 
mean transition time. The thermal fluctuation can be 
neglected only when the transitions take place in a time 
much larger than the observational time. 

In all the previous situations the quantum theory ap- 
plies. Because the density matrix automatically includes 
the mean over an ensemble, there is no clear distinction 
between the dissipative dynamics and the dynamics in- 
duced by the fluctuations. QSD, by unraveling the dif- 
ferent quantum trajectories, helps one to understand the 
role played by the statistical mean. 

VI. SUMMARY 

We have shown that QSD leads to quantum trajecto- 
ries which exhibit bistability for the driven damped an- 
harmonic oscillator and that the quantum steady state 
result is recovered through random switching between the 
two equilibrium points due to quantum fluctuations. The 
fluctuations also introduce a characteristic time scale: 
the mean transition time. 

We have also shown that it is still possible to observe 
bistability dynamically in this quantum system by in- 
troducing a measurement delay t m . An experiment will 
show hysteresis only if the transition time between the 
(classical) equilibrium points, r, is much larger than all 
other relaxation times involved (approximately 1/k). If 
the transition time is not that large, it is not possible 
to observe any hysteresis effects and the quantum steady 
state solution is expected. Furthermore, provided the 
transition time is much larger than the other times, hys- 
teresis can be seen only if a measurement delay is such 
that t m <C t. Thus, within this interpretation the clas- 
sical result is valid only if the quantum fluctuations are 
so small that they induce a transition time very large 
compared to the period of observation. 

Our results confirm the statistical description given by 
Drummond and Walls. The results of Bortman and Ron 
have to be examined more carefully. Strictly speaking, 
even for low excitation number, the bistability is not de- 
stroyed, it is just hidden by the quantum fluctuations. 
But one cannot neglect the effect of the fluctuations upon 
the two stable states since it is exactly these fluctuations 
which prevent an experimental observation of hysteresis 
in the steady state. 

Finally, if QSD is used to describe the dynamical be- 
haviour of a single quantum system our investigations 
might be of particular relevance for experiments on opti- 
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cal bistability or for the recently performed experiments 
of Gabrielse et al. 0] in which hysteresis was observed. 
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FIG. 1. Steady state excitation number \a\ 2 (classical re- 
sult - dots) and (o'a) (quantum result - line) versus the de- 
tuning Au>. The parameters are: the damping k = 1.5, the 
driving (3 — — 7.0 (choosing f3 real is equivalent to setting the 
origin for the phase of (a), introducing a complex /3 simply 
produces a rotation in phase space) and the anharmonicity 
X = 0.05. 

FIG. 2. Representation of the domain of hysteresis. The 
three conditions given by Eqn.(^) in the text describe the 
border of this domain. Inside the bounded region the system 
possesses three solutions, two stable and one unstable, this is 
the bistable domain. Outside, the system has only one solu- 
tion, always stable. Notice that the bistable region is entirely 
defined by only two parameters x = A J +X and y — l'i ■ 



FIG. 3. Evolution of a quantum trajectory in phase space. 
The parameters are: k, = 1.5, /3 = —7.0, x = 0-05 and the 
detuning Acu = —5.0. The initial state is chosen to be a 
coherent state centered at ((Q) = 7, (P) = 14). The first stage 
of the evolution is the decay towards a local minimum, in a 
time of the order of 1/k. Then the system fluctuates around 
the equilibrium point for an amount of time given, in mean, by 
the transition time. Finally, a big enough fluctuation occurs 
to project the system into the other basin of attraction where 
the system remains for a very long time. 

FIG. 4. Representation in time of A) the position (Q) (full 
line) and its variance AQ 2 (dotted line) and B) the momen- 
tum (P) (full line) and its variance AP 2 (dotted line) at 
the particular instant of the transition (approximatively at 
t = 596.3 in this example). Same parameters as figure [| 
Notice the derealization in space of the wave packet at the 
transition. Before and after the transition, the variances AQ 2 
and AP 2 are small (compared to the distance between equi- 
librium points) showing localized states. 

FIG. 5. Evolution in time of the mean position (Q) (lower 
curve) and momentum (P) (upper curve). Parameters same 
as figure ||. The mean is taken over 100 realizations. The time 
scale of the decay is much larger than the dissipative time of 
1/k = 0.67. 

FIG. 6. Simulation of an ideal (single) experiment accord- 
ing to the QSD model. Parameters are the same as in fig- 
ure ^. The detuning step is 0.1 and the measurement time is 
t m = 50. The dotted line represents the classical steady state 
excitation number. 
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